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ABSTRACT 

Hamilton et al. (1991) proposed a simple formula relating the nonlinear autocorre- 
lation function of the mass distribution to the primordial spectrum of density fluc- 
tuations for gravitational clustering in an £1 = 1 universe. We present high resolu- 
tion N-body simulations which show this formula to work well for scale-free spectra 
P(k) oc k n when n ~ 0, but not when n < — 1. We show that a modified version of 
the formula can work well provided its form is allowed to depend on n. Our modi- 
fied formula is easy to apply and is a good fit to our N-body simulations which have 
< n < — 2. It can also be applied to non-power law initial spectra such as that of 
the Cold Dark Matter model by using the local spectral index at the current nonlinear 
scale as the effective value of n at any given redshift. We give analytic expressions 
both for the nonlinear correlation function and for the nonlinear power spectrum. 
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INTRODUCTION 



According to the theory of gravitational instability, the clus- 
tering of matter in the universe is determined by the power 
spectrum of primordial density fluctuations. Since present 
cosmic structures are highly nonlinear on small scales, the 
relation between primordial and present spectra results from 
the interplay of complex processes. It is nevertheless highly 
desirable to find a simple formula which can relate the two 
at least approximately. Such a formula not only gives the 
two-point clustering properties of the evolved density field 
for any given initial model, but may also enable us to recon- 
struct the primordial density spectrum from the observed 
correlations of galaxies. 

It is generally assumed that in an Einstein-de Sitter 
universe, structure will grow self-similarly in time if the pri- 
mordial spectrum has a power law form: 



P(k) = Ah 71 



(1) 



where A is a constant and 4 > n > —3. For a given spectral 
index n, the correlation functions at different epochs are 
then related by a simple scaling relation (see Efstathiou et al. 
1988, hereafter EFWD). Realistic power spectra are not pure 
power laws, but we will show that their evolution between 
the linear and highly nonlinear regimes is similar to that of 
a power law spectrum with an effective index n e s where we 
define n e e by 



where ro is the radius of the top-hat window in which the 
rms mass fluctuation is unity. 

Hamilton et al. (1991, HKLM) made an ansatz for the 
relation between the linear, and the evolved (or nonlinear) 
average two-point correlation functions £l and £e: 

£e(P) = F[£l(Rq)] , Ro = [1 + &(R)] 1/3 R. (3a) 

Here f(r) = (3/r 3 ) J y 2 £,{y)dy, and F is a universal func- 
tion assumed to be independent of the initial spectrum. This 
ansatz relates £e at a given length scale to £l at a larger 
scale, thus taking account of the fact that nonlinear den- 
sity fluctuations shrink in comoving coordinates as evolu- 
tion proceeds. A remarkable aspect of this ansatz is that 
the same function F appeared to work for all epochs and 
initial spectra (see also Nityananda & Padmanabhan 1994). 
HKLM showed that their ansatz is in good agreement with 
the N-body simulations of EFWD, and used these data to 
determine the functional form of F. Motivated by this work, 
Peacock & Dodds (1994, PD) studied a similar ansatz for 
the dimensionless power spectrum: 

" 1/3 fc, 



A E (fc) = #[A L (*o)] , k = [1 + A E (fc)r 



(36) 



Icff 



dlnP(fco) 



din fci 



(2) 



where Ae(A;) = 4-7rfc Jfc(fc) (with Pe being the evolved 
power spectrum), Al(A:o) = 47rfcoP(fco), and $ is a universal 
function. The relations in equations (3) are independent of 
the form of the initial spectrum P, but the authors noted 
that for n = —2 their formulae were not in good agreement 
with the simulation data. Such formulae are extremely use- 
ful, so it is important to check their validity for a wide range 
of initial spectra by using numerical simulations with better 
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resolution than the old EDFW models. 

This paper provides such a check. We find that equa- 
tions (3a) and (3b) both work remarkably well for power law 
spectra with n ~ 0, but that they are seriously in error for 
n < — 1. This deficiency can be rectified by making the func- 
tional forms of F and $ depend on n. We provide improved 
formulae which include a simple n— dependence and show 
how they can be applied to non-power law models such as 
CDM. Section 2 presents results from several N-body simu- 
lations and compares them with the formulae of HKLM and 
PD. We present our improved model in Section 3 and give 
our conclusions in Section 4. 



spectra. 

The dotted curves in Figure 2 show Ae(&) as a function 
of Al(A;o) for the same simulations as in Figure 1. The self- 
similarity of clustering is again verified for the power law 
models, though the results for n = —2 have some scatter. 
Again the relation between Ae(&) and Al(£:o) depends on n. 
The long dashed curves in each panel show the fit proposed 
by PD. As for the HKLM formula, the PD formula works 
well only for n ~ 0. For n < — 1 and for the CDM spectrum, 
there are substantial deviations from the simulation data. 
This is not surprising because, as pointed out by PD, their 
formula reproduces the HKLM formula very accurately for 
fl = 1, and so reproduces its shortcomings. 



2 NONLINEAR EVOLUTION AND THE 
UNIVERSAL SCALING ANSATZ 

In this paper, we give simulation results for six different 
spectra. One is the standard CDM spectrum. The other 
five are power-law spectra with n = 0, —0.5, — 1, —1.5 and 
—2. All simulations assume an Einstein-de Sitter universe. 

These simulations were performed using high resolu- 
tion particle-particle/particle-mcsh (P 3 M) codes. For mod- 
els with > n > —1.5, the code is the same as in EFWD, 
but was run with more particles (100 3 , compared to 32 3 ) 
and higher force resolution (L/2500, compared to L/200, 
where L is the side of the computational box). The power 
spectrum is normalized as described in EFWD. 

The n = —2 run was performed by E. Bertschinger us- 
ing an adaptive P 3 M code. It followed 128 3 particles with a 
force resolution of L/2560. The CDM simulation, performed 
by Gelb & Bertschinger (1994), is also a P 3 M simulation. It 
has 144 3 particles, Q = 1, H = 50km/s/Mpc, L = 100 Mpc 
and a force resolution of 65kpc. It is normalized so that as 
(the linear rms mass fluctuation in a sphere of radius 16Mpc) 
is unity when the expansion factor a = 1. For the n = — 2 
and CDM simulations we have checked our results against 
other independent simulations — for n = —2 against an en- 
semble of three N = 256 3 PM simulations, and for CDM 
against an N = 100 3 P 3 M simulation in a 50 Mpc box. In 
all cases the results agree to better than about 20% in the 
critical transition between linear and nonlinear regimes (see 
below) . 

The dotted curves in Figure 1 show £e(R) as a function 
of £l(Ro)- In each case we show results for four different 
expansion factors as detailed in the figure caption. For each 
n, the curves for different a are close to each other, thus 
demonstrating that clustering is self-similar in the power 
law models (see Jain & Bertschinger 1995 for a detailed ex- 
amination of self-similar scaling for —3 < n < —1). The 
flattening of the curves at high amplitudes, which is severer 
for earlier output times, reflects the small scale resolution 
limit. These results show that the log^E(i?)-log^L(i?o) re- 
lation depends on n and is a good fit to the HKLM formula 
(the long-dashed curve in each panel) only for n ~ 0. It 
misses the simulation results by factors of up to 3 and 10 
around £l ~ 1 for n = —1.5 and —2 respectively. A sim- 
ilar failure is also seen in the CDM case where results for 
different a (or erg) differ substantially from each other. At 
a — 0.2 the discrepancy around £l ~ 1 is almost a factor of 
10. The effective power indices n e ff are —0.7, —1.3 and —2 
for a — 1, 0.5 and 0.2 respectively. The dependence on n e s 
is strikingly similar to that for the corresponding power-law 



3 AN IMPROVED MODEL 

Since for given n the correlation functions and dimensionless 
power spectra obey scaling relations of the form (3), we can 
write both in the more general form y = f n (x). From Fig- 
ures 1 and 2 we see that although the relation between the 
evolved and unevolved quantities differs for different spectra, 
the shapes of the various curves are similar. This suggests 
that the n-dependence can be scaled away by a simple shift 
in the log-log plane. From linear theory and from the stable 
clustering hypothesis respectively we know that /„ has the 
following asymptotic behaviour: f n (x) — > and 
f n (x) oc x :i ^ 2 as x — » co. Thus the shift can be made only 
in the direction y = x. This suggests the ansatz: 



faQR) 

B e (n) 

and 

A E (fc) 
B A {n) 



= F 



= <E> 



ClORo) 



B ( (n) 



A L (fc ) 



B A (n) 



(4a) 



(4b) 



where B^(n), B A (n) are constants depending on n, and F 
and $ remain independent of n. Since F and $ are non- 
linear functions, B^{n) and Ba(it,) are generally not equal, 
although A and £ are proportional to each other. The sim- 
ulation data can be fitted reasonably well by the following 
simple empirical formula: 



Be(n)=(—J ; 



(5a) 



(5b) 



As n approaches —3, both -Bj(n) and -Ba(^) tend to zero 
and our approach breaks down. This reflects the breakdown 
of hierarchical clustering as n — > —3. 

In the first panel of Figure 1, the solid curve shows our 
best fit to the simulation results for n = 0. Our result is 



F(x) 



x + 0.45a; 2 - 0.02a: 5 + 0.05a: 6 



(6a) 



1 + 0.02a; 3 + 0.003a; 9 / 2 

The inverse of this can be approximated to better than 4% 
over the range of interest by 

1/3 



f(v) = v 



1 + O.OSet; 1 - 93 + O.OOOlt; 3 
1 + 1.75j/ - 0.0015i/ 3 - 63 + 0.028y 4 



(6b) 



The solid curves in the other panels of Figure 1 are obtained 
from equation (4a), with F given by (6a) and B^(n) given 
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by (5a). For the CDM case, we have used the effective power 
indices given above. It is remarkable that this simple for- 
mula works accurately in all cases. For spectra with n < — 1 
and for the CDM spectrum, it makes a substantial improve- 
ment over the formula given by HKLM. For all six spectra, 
the N-body points (the mean of the 4 time outputs in the 
case of n = — 2) agree with our fit to better than 20% for 
£e ~ 1 — 100, provided it is within the range of numeri- 
cal validity of the simulation for the given time. Since the 
scatter between different simulations of a given spectrum is 
typically 15%, this fit is sufficiently accurate. 

The solid curve in the first panel of Figure 2 is given by 
the fitting formula 

_ , f l+O.Gx + x 2 -0.2x 3 -1.5x 7/2 +x 4 \ 1/2 ,„ . 
*<*> = X { 1 + 0.0037,3 ) ■ ^ 

The inverse of this function is fitted to an accuracy of better 
than 2% over the range of interest by 

t(y)=,f l + 0^- 86 + 0.000345^ V /3 
\y> y ^i + i. 582y + o.05j/ 3 +0.072y 4 y v ; 

The solid curves in other panels are obtained from equation 
(4b), with $ given by (7a) and B&{ri) given by (5b). For 
the CDM spectrum, we have again used the effective power 
index n e g. We have transformed the AE(fc)-AL(fco) rela- 
tions represented by the solid curves to the corresponding 
£e{R)-£l(Ro) relations. In all cases the transforms agree 
with the £e{R)-£l{Ro) relations shown by the solid curves 
in Fig.l. 

The figure shows that our model also works for the di- 
mensionless power spectrum. There are small differences in 
the shapes of the N-body curves for different spectra — as n 
becomes more negative, the AE(fc)-AL(fco) curve rises more 
sharply in the regime where Al(£;o) ~ 1. For CDM the effect 
is pronounced at high redshifts. This suggests that one can 
do marginally better by introducing an n— dependence in the 
shape of the curve as well. There is also a hint that for CDM- 
like spectra the curve steepens more than the corresponding 
power law case. However our fit is adequate given the scat- 
ter in the simulation data, as the discrepancy between the 
fit and the simulations rarely exceeds the level of the scatter 
in the simulation data, which is typically 15 — 20%. This 
agreement can be compared with the discrepancies of the 
PD fit, which is almost a factor of 10 below the mean of the 
simulation data for n = —2 and for CDM at a — 0.2 around 
A E (fc) ~ 10. 

We note that the lowest order nonlinear term in the fit- 
ting formulae (6a) and (7a) is of order x 2 . However, this 
term is significant only in the range 1 < £l,Al < 4 - at 
smaller amplitudes linear theory is sufficiently accurate, and 
at larger amplitudes higher order terms dominate the fitting 
formulae. Whether second order perturbation theory agrees 
with our formulae within the range of dominance of the x 2 
term can be ascertained by comparing the coefficient of this 
term with the predictions of perturbation theory. The out- 
come is dependent on the spectrum, because perturbation 
theory does not treat different spectra in the same way as 
the scaling ansatz. For non-power law spectra the outcome 
also depends on the time at which the comparison is made; 
for example, for the CDM spectrum the agreement is very 
good at high redshift, but for z £ 1 perturbation theory 



significantly over-estimates the nonlinear spectrum. Since 
our formulae fit the simulation results, these features can 
be inferred from previous work in which the second order 
power spectrum is compared with N-body simulations (Jain 
& Bertschinger (1994), and references therein). 



4 DISCUSSION 

For a given linear power spectrum P(fco), the effective power 
index n e s can be obtained from equation (2). The linear di- 
mensionless power spectrum AL(fco) and the linear average 
correlation function £l(-Ro) can be calculated directly from 
P(fco). The evolved dimensionless spectrum Ae(/c) and the 
evolved average correlation function £e{R) can then be cal- 
culated from equations (4b) and (4a), respectively, with F 
given by equation (6a), $ given by equation (7a), and B(n) 
given by equations (5). To obtain the linear spectrum from 
the evolved spectrum (or correlation function), one can use 
equations (6b) and (7b). In this case, an iterative procedure 
may have to be used, because the effective power index is 
unknown a priori. As noted by PD the effects of redshift 
space distortions and biasing also need to be considered in 
such reconstructions. For other applications, such as using 
the nonlinear power spectrum to compute the statistics of 
photon trajectories, our results can be used directly (Seljak 
1995). It would clearly be useful to extend these results to 
cosmological models with SI < 1 or which involve a hot dark 
matter component (see PD). 
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FIGURE CAPTIONS 

Figure 1. The evolved average mass correlation function 
£e(-R) as a function of the linear average mass correlation 
function £l(#o)- Note that these two functions are calcu- 
lated at two different scales, as discussed in the text. The 
dotted curves show the results derived from N-body simula- 
tions of different spectra. Results are shown for four different 
expansion factors a; a curve that flattens earlier corresponds 
to a lower value of a. The expansion factors are a = (6.10, 
14.91, 36.86, 90.88) for n = 0, (4.52, 9.53, 20.22, 42.82) for 
n = -0.5, (3.34, 6.08, 11.08, 20.20) for n = -1, (2.47, 3.86, 
6.07, 9.52) for n = -1.5, (1.14, 1.62, 2.29, 2.72) for n = -2, 
and (0.2, 0.5, 1.0) for CDM. Long-dashed curves show the 
fitting formula of HKLM. Solid curves show the results of 
our improved formula (see text). The two light dotted lines 
show the expected behaviour of the relation in the linear 
and stable clustering (highly nonlinear) regimes. 

Figure 2. The evolved dimensionless power spectrum 
Ae(&) as a function of the linear dimensionless power spec- 
trum Al(&o)- Note that these two functions are calculated 
at two different scales, as discussed in the text. The dotted 
curves show the results derived from the same N-body sim- 
ulations as described in Figure 1. Long-dashed curves show 
the fitting formula of PD. Solid curves show the results of 
our improved formula (see text). The two light dotted lines 
show the expected behaviour of the relation in the linear 
and stable clustering (highly nonlinear) regimes. 
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